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Abstract. We calculate upper and lower bounds on the modulus of the pion electromagnetic form factor on 
the unitarity cut below the ljtv inelastic threshold, using as input the phase in the elastic region known via 
the Fermi- Watson theorem from the mr P-wave phase shift, and a suitably weighted integral of the modulus 
squared above the inelastic threshold. The normalization at t = 0, the pion charge radius and experimental 
values at spacelike momenta are used as additional input information. The bounds are model independent, 
in the sense that they do not rely on specific parametrizations and do not require assumptions on the phase 
of the form factor above the inelastic threshold. The results provide nontrivial consistency checks on the 
recent experimental data on the modulus available below the ljtt threshold from e + e~ annihilation and 
r-decay experiments. In particular, at low energies the calculated bounds offer a more precise description 
of the modulus than the experimental data. 
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1 Introduction 

The pion electromagnetic form factor F(t) is a fundamen- 
tal observable that has been used as a probe of QCD 
in several regimes, including Chiral Perturbation Theory 
(ChPT) and lattice simulations at low energies and per- 
turbativc QCD at asymptotic spacclikc momenta. On the 
other hand, even before the rise of QCD as the modern 
theory of strong interactions, form factors were analyzed 
based on the general principles of analyticity and unitar- 
ity. 

The experimental information on the pion form fac- 
tor improved significantly in recent years. In particular, 
the modulus |P(t)| on the unitarity cut t > £ + , where 
t + = 4M^, was measured by several e + e~ — > tt + tt~ ex- 
periments 

0, 0, H II Hi and, invoking isospin symme- 
try, also from hadronic r-decays @. In the elastic region 
t + < t < ii n where t- m = (M u + M^o) 2 is the first sig- 
nificant inelastic threshold, the phase is known via the 
Fermi- Watson theorem from the phase-shift of the P-wave 
7T7T elastic scattering amplitude, determined recently with 
great accuracy from Roy equations and scattering data 
0, IE H: EH EH • While an obvious correlation between the 
phase and modulus is implied by explicit formulae such 
as Breit-Wigner parametrizations, independent sources of 
correlation are required to test these measurements in a 
model independent way. 

In this work we present upper and lower bounds on 
the modulus of the form factor in the clastic region, re- 



sulting from phase information on the unitarity cut in the 
clastic region. Additional inputs coming from the normal- 
ization at t = 0, the knowledge of the pion charge radius 
[T3 . HH] and spacelike measurements [lj, EH can offer fur- 
ther constraints. The appropriate tools for this investiga- 
tion are analyticity and unitarity, which correlate the var- 
ious pieces of the input. Recently developed mathematical 
techniques [13, 13 ) which were used to constrain the 
pion form factor near the origin [l9| and on the spacelike 
axis H3, can be readily extended to the elastic region of 
the unitarity cut. 

The method uses as first input the knowledge of the 
phase 



Aig[F(t + ie)] = 5{(t), t+<t<t h 



(1) 



where 8\ (t) is the phase shift of the P-wave of irir elastic 
scattering. 

The second piece of input comes from the knowledge 
of \F(t)\ above i; n . The modulus has been measured up to 
\ft = 3 GcV by BaBar Collaboration Combined with 
the asymptotic decrease \F(t) \ ~ 1/t at large t > 0, which 
can be inferred from the perturbative QCD prediction at 
large space like momenta and mathematical theorems on 
the asymptotic behaviour of analytic functions [2l|, this 
information can be used to calculate an integral defined 



1 



dtp(t)\F(t)f 



I. 



(2) 
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Here p(t) is a positive-definite weight, for which the inte- 
gral converges and an accurate evaluation of / is possible. 
The optimal procedure is to vary p(t) over a sufficiently 
large admissible class and take the best result. A suitable 
class is given by the simple expression 



Pit) 



(3) 



where Q 2 > and b, c taken in the range b < c < b + 2. 
The dimension of / can be inferred from Eqns.© and ©. 
As we shall discuss, in the choice of Q 2 , b, c a compromise 
should be achieved between the objective of deriving suffi- 
ciently strong bounds and the need of a precise evaluation 
of the integral © with the present knowledge of |.F(t)|. 

Additional information inside the analyticity domain 
can also be implemented exactly. In practice we shall use 
the normalization and the first derivative at t = 0: 



F(0) = 1, 



F'(0) = l(rl 



(4) 



with the charge radius (r 2 ) varied within reasonable limits 
[HI, [l3| . We shall further implement the value of the form 
factor at a spacelike momentum t s < 0, 



F(t s ) = F s ±e s 



(5) 



known from the most precise experiments [14 . If 5| . 

The relations © , © , © and © define a class of 
functions F(t) real analytic in the i-plane cut for t > t+. 
In Refs. [f9|, [2(j, bounds on the higher Taylor coefficients 
at t = and on the values of F(t) on the spacelikc axis, 
for functions F(t) belonging to this class were derived. We 
extend now the formalism and derive rigorous upper and 
lower bounds on the modulus \F(t)\ in the region t+ < 
t < tj n , by solving the corresponding extremal problem on 
the same functional class. These bounds will be compared 
with recent experimental measurements of the modulus 
from e + e~ — > tt + tt~ [H, H, H 0, Q and r-decays @ below 
the W7T threshold. 

We emphasize that the goal of present work is to check 
(and confirm) the consistency of various data sets on the 
modulus with the bounds derived from the input described 
above. Therefore, we have not used as input in the formal- 
ism the values of the modulus below the inelastic thresh- 
old provided by the recent experiments. Adding to the 
input one or more of these measured values would lead 
to a considerable strengthening of the bounds. However, 
as the different sets of data are not perfectly consistent 
amongst themselves, the choice of the input from one set 
would introduce a certain bias, which we prefer to avoid 
in this first analysis. It may also be noted that one could 
use the techniques applied here to improve, for instance, 
the extracted value of the pion radius, using the precise 
values of the recently measured modulus. These problems 
are of interest and will be considered in a future work. 

The scheme of the paper is as follows: in Sec. [21 we give 
a brief review of the mathematical formalism. In Sec. [3l 
we describe the information used as input, discussing also 
the inclusion of the main isospin breaking effect, produced 



by the p — uj interference, which affects the form factor ex- 
tracted from e + e~ annihilation. In Sec. [4l we present the 
results of our computations, including a detailed analysis 
of the uncertainties of the input, and make a comparison 
with the most recent experimental measurements of the 
modulus from e + e~ — > 7r + 7r~ " 0, i 1 1, d and r-decays 
6]. Finally in Sec. [5j we provide a discussion and a sum- 
mary of our results. 



2 Method 

We exploit © by introducing the Omnes function defined 

by 



0(t) = exp 



t 



dt 



, 5(t<) 



t'(t> -t) ' 



(6) 



where S(t) = S\(t) for t < t- m , and is an arbitrary function, 
sufficiently smooth (i.e., Lipschitz continuous) for t > t- m . 
The crucial remark is that the function h(t), defined 

by 

F(t)=0(t)h(t), (7) 

is real for t < tj n , i.e. it is analytic in the t-plane cut 
only along t > t; n . Moreover, from © it follows that h(t) 
satisfies the condition 



1 



dtp(t)\0(t)\ 2 \h(t)\ 2 =L 



(8) 



'in 



In order to exploit this property, we cast the problem into 
a canonical form by performing the conformal transforma- 
tion 



tin y^in ^ 



Un + \Ain — t 



(9) 



which maps the complex t-plane cut for t > t; n onto the 
unit disk \z\ < 1 in the z-plane defined by z = z(t). The 
upper (lower) lip of the branch-cut [t- m , oo] is mapped onto 
the upper (lower) half of the unit circle in the complex z- 
plane, the real line [— oo,0] to [—1,0] and [0,<i n ] to [0,1]. 

The next step is to introduce two outer functions, 
i.e. functions analytic and without zeros in the unit disk 
\z\ < 1, defined in terms of their modulus on the bound- 
ary, related to y/ p(t) \dt/dz(t)\ and |0(f)|, respectively 
[la , In particular, for weight functions with the ex- 
pression ©, the first outer function w(z) can be written 
in an analytic closed form in the z- variable as [lTj 

i+6-c (1-^) 1/2 [l + z{-Q 2 )Y 

(l + 2)3/2-c+& (! _ Z Z{-Q 2 )Y 



W(z) = (2v^n) 



(10) 

where z(— Q 2 ) is computed using Eq. ©. 

For the second outer function, denoted as w(z), we 
shall use an integral representation in terms of its modulus 
on the cut t > t- m , which can be written as [TH, [TtJ 



u>(z) = exp 



\n\0(t')\dt' 



tin Vr^Mf -*(*)) 



(11) 
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where t(z) is the inverse of z = z(t), for z(t) defined in 
([9]). We mention that the product w(z)co(z) defines in 
fact a single outer function with modulus accounting for 
all the factors in ([5]). The separation into two functions 
was suitable in our case since one of them, w(z), could be 
written in the simple analytic form (TTU)) . 
Further, if we define a function g(z) by 



g(z)=w(z)uj(z)F(t(z))[0(t(z))] 



(12) 



it is easy to see that the integral (0) can be written in 
terms of |.g(z)| 2 for z on the boundary (z — e ) as 



l - de\g{e i6 )\ 2 = I. 
■~ Jo 



(13) 



The £ 2 -norm condition (fT3|) leads to what is known as the 
Meiman interpolation problem, which consists in finding 
the most general rigorous correlations between the val- 
ues of the function and its derivatives inside the unit disk 
\z\ < 1. These conditions can be written either as a posi- 
tivity of a determinant and its minors, or as an inequality 
for an explicit convex quadratic expression of the input 
quantities (for a proof and older references see [17| ). For 
instance, if we denote go = <?(0) and g\ = <?'(0) the value 
of the function g(z) and the first derivative at z = 0, 
and by g(z a ) and g(zb) the real values at two real points 
\z a \ < 1 and \zb\ < 1, then (fl3|) implies the determinantal 
inequality 



I-9o~ 9i 9(z a ) g(z b ) 



g{z a ) 



Zj (ZgZ b ) 2 

1 - zl 1 - z a z b 
{z a z b ) 2 (z b ) 4 



1 - z a z b I - zl 



> o, 



(14) 



where 



We recall that \0(t)\ is obtained from (jSJ) via the Principal 
Value (PV) Cauchy integral 



g(z a ) = g(z a )-g -giz a , g(z b ) = g{z b )~ g - g x z b . (15) 

Furthermore, all the minors of the determinant above should 
be nonncgative. 

In the present case we note that ^ implies that the 
origin t = of the i-plane is mapped onto the origin z = 
of the z-plane. Then the values go and g\ are related in 
a straightforward way to F(Q) = 1 and the charge radius 
(r 2 .). Also, we shall take z a = z(t a ) < 0, with t a = t s a 
spacelikc point included as input, and < z b = z(t b ) < 
1, where t b is an arbitrary point in the range (< + , tin). 
Then using as input F(t a ) = F(t s ) from ([5]), i.e. a known 
value for 

g{z a ) =w(z a )u)(za)F(t a )[0(t a )]- 1 , z a = ~z{t a ) 1 (16) 



the inequality (|14[) leads to upper and lower bounds for the 
real values g(z b ). From these we obtain upper and lower 
bounds on the modulus |F(if,)|, written using ([12"]) as 



\nt b )\ = \o(t b )\ 



g{z{h)) 



w{z{t b )) u(z{t b )) ' 



(17) 



\(D(t)\ = cxp ( -PV [ 

\7T Ji 



dt' 

4M 3 t'(t' - t) 



(18) 



It is important to emphasize that, as proven in [16l . ll7| . the 
bounds are independent of the phase 8(t) for t > t\ „ en ter- 
ing the Omnes function ([5]). Also, it is easy to show |16l.[l~7| 
that, at fixed p(t), the bounds depend in a monotonous 
way on the value of the quantity /, becoming stronger 
when this value is decreased. 



3 Input 

Recall that the first inelastic threshold t- m is due to the 
opening of the ujt: channel, ^Jt^ = 0.917 GeV. Below i in 
we use as phenomcnological input the phase shift 5{ (t) de- 
termined recently with high precision from Roy equations 
satisfied by the tttt elastic amplitude @, H, H, E3, E3 ■ In 
particular, the P-wave phase-shift is parametrized in fiol ] 



cot5l(t) = £ s (M* p - 



■t) 



( 2Mg 
\AP p Vt 



Bo + Bi 



Vt-Vh 



where k = y/t/4 - M% and y/to = 
0.011, B 1 = 0.19 ± 0.05 and M p 



(19) 

1.05 GeV, B = 1.043± 
= 773.6 ± 0.9 McV (we 
have taken the parameters from the so-called CDF solu- 
tion [13). 

The phase 5\(t) obtained from (fH?)) with the central 
values of the parameters is very close to the solution of 
the Roy equations presented recently in [ll] ]: the differ- 
ences between the two phases are less than a small frac- 
tion of a degree at low energies and reach one degree only 
at the upper end t- m of the elastic range. In fact, these 
differences are comparable with the uncertainty resulting 
from the errors of the parameters quoted below (|19p . In 
our analysis we have used for the central value of the phase 
the expression p^|) and estimated the uncertainty either 
from the errors of the parameters quoted above or as the 
difference between the phase (|19p and the phase obtained 
in 
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The phase shift given above was obtained assuming 
isospin symmetry, with = . As discussed in [22] , 
an important isospin breaking effect, which manifests itself 
in e + e~ annihilation, is produced by p — uj interference. 
We draw attention also on the extensive analyses of the 
data from e + e~ annihilation and r-decay, in particular of 
the isospin breaking effects, performed recently in [23l.l2~H. 
while the comparison of experiment and theory in e + e~ — > 
hadrons was discussed in detail in [25|, [2(| , in connection 
with the hadronic contribution to the anomalous magnetic 
moment of the muon (g — 2). 

In our work we shall model the isospin breaking pro- 
duced by p — uj interference effect by adding to the phase 
8(t) entering the Omnes function © the correction A8(t) = 



4 
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Table 1. Values of the integral / defined in (J2]| corresponding 
to the weight functions adopted in our analysis. Note that 



all energies are 


in GeV. 






b 


c 


Q 2 


I 





1 





0.578 ± 0.022 


1/2 


1 


3 


0.246 ± 0.011 



Arg [F p _ u (t)] , where 
,(*)= (l + e- 



F n 



C = (M W -V2^) 2 , (20) 



with M u = 0.7826 GcV, r u = 0.0085 GeV, and the small 
parameter e = 1.9 X 10~ 3 . 

Alternatively, the effect can be accounted by multi- 
plying the bounds derived in the isospin limit with the 
modulus of the correction factor Fp^^it) given in (l2"0l) . 
The two procedures give practically the same results. 

Above <; n = (0.917 GeV) 2 we use in © a continu- 
ous function 8{i), which approaches asymptotically ir. As 
shown in [l7j . if this function is Lipschitz continuous, the 
dependence on S(t) of the functions 0(t) and uj(z), de- 
fined in ([6]) and pip , respectively, exactly compensate 
each other, leading to results fully independent of the un- 
known phase in the inelastic region. This is one of the im- 
portant strengths of the method advocated in this work. 

For the calculation of the integral defined in we 
have used the BaBar data [l[ from <; n up to \fi = 3 GeV, 
continued with a constant value for the modulus in the 
range 3 GeV < \fi < 20 GeV, smoothly connected with 
a 1/i decrease above 20 GeV. As discussed in [2(J these 
assumptions are expected to overestimate the value of /. 
In view of the monotony property mentioned at the end of 
Sec. [21 this provides a guarantee of robust and conservative 
bounds. 

We have worked in our analysis with weights of the 
type (J3J) ■ As shown in [2C$, for wcigths decreasing as l/\/t 
or faster, the contribution of the asymptotic regime \fi > 
20 GeV is negligible. We have tested a large class of ex- 
pressions, and finally settle down to two suitable weights: 
p(t) = l/t and p{t) =y/i/(t + 3). The values of / for the 
corresponding choices of the parameters b, c and Q 2 are 
given in Table [TJ where the uncertainties are due to the 
BaBar experimental errors. 

As discussed in Sec.[I] we use as input also some values 
inside the analyticity domain, as specified in (g]) and ([5]). 
For the charge radius we have adopted the range given in 

dm 

(r 2 ) =0.43 ±0.01 fm 2 , (21) 
while for the spacelike data we have taken the values 



F(-1.60 GcV 2 ) = 0.243 ± 0.012±g;gJ| , 



0.167 ±0.010 



+0.013 
0.007 ) 



F(-2.45 GeV 2 
given by the most precise recent experiments [1 



(22) 



30 



20 



10 



— no spacelike 

— spacelike 



~i 1 r 



p(t) = 1/t 




0.5 0.6 0.7 

t" 2 [GeV] 



Fig. 1. Upper and lower bounds on |F(t)| 2 below 0.917 GeV for 
central values of the input, obtained with the weight p(t) — l/t. 
Dashed lines: bounds obtained without spacelike input; solid 
black lines: bounds obtained with spacelike input. 
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Fig. 2. Comparison of the two methods of including the un- 
certainties of the input discussed in the text. 
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Fig. 3. Comparison of the two methods of accounting for the 
errors on the phase shift <5i(t) discussed in the text. 



4 Results 

In Fig. [1] we show the upper and lower bounds on |-F(£)| 2 as 
function of the centre of mass energy y/t, on the unitarity 
cut up to the ujit threshold \ft\[\ = 0.917 GeV, derived 
with the formalism presented in Sec. [H The curves are 
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obtained with the central values of all the input quantities 
described in Sec.[3l in the isospin symmetry limit, using in 
(J2J) the weight p(t) = 1/t. The value / of the integral for 
this weight is given in Table [TJ At each energy, the allowed 
values of \F (t)\ 2 are situated between the two solid black 
curves. 

The figure illustrates also the effect of including a space 
like datum: the dashed red lines are the weaker bounds 
obtained without the spacelike information. As discussed 
in [20j |. in order to improve the bounds on F(t) at large t 
on the spacelike axis, it was convenient to use as input the 
value of F(t) at a larger \t\, i.e. the second value given in 
(22). For the bounds on the modulus on the timelike axis 
the effect is not so clear-cut and depends on the energy, 
but in general it turns out that the first value in (f2"2")l is 
useful especially for improving the upper bounds, while 
the second is more efficient for the lower bounds. There- 
fore, the black solid lines in Fig. [T] arc obtained using as 
input the value of F(-1.60 GeV 2 ) for the upper bounds 
and the value F(— 2.45 GeV 2 ) for the lower bounds, and 
this choice will be made throughout the subsequent anal- 
ysis. 

The uncertainties of the input are expected to lead to 
weaker bounds and to enlarge the allowed range for the 
modulus at each energy. In order to investigate this effect 
we have used two methods: 

— Method 1 - vary separately each input inside its error 
interval with the others kept fixed at their central val- 
ues, and then combine the corresponding errors on the 
bounds in quadrature. 

— Method 2 - vary simultaneously all the input variables 
inside their allowed intervals set by the quoted errors, 
and take the most conservative results, namely the 
largest upper bounds and the smallest lower bounds 
yielded by a compatible set of input quantities inside 
these intervals. 




t 1 ' 2 [GeV] 



Fig. 4. Bounds on |.F(i)| 2 below t- m obtained with the weight 
p(t) = 1/t in the isospin limit. The central grey band is defined 
by the upper and the lower bounds obtained with the central 
values of the input quantities, the extended cyan bands are 
obtained by varying all the inputs within the error intervals. 
Also shown are the Belle data [al. 
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• Belle | 1 




P(t) 


= 1/t 


i 


i.i, 



1.3 0.35 0.4 0.45 0.5 

t 1 ' 2 [GeV] 



Fig. 5. Low-energy profile of the bounds obtained with the 
weight p(t) = 1/t in the isospin limit. For details see caption 
to Fig. H 

We have found that these two methods give similar 
results, as illustrated in Fig. [3J where the bounds are ob- 
tained with the weight p(t) = 1/t in the isospin limit. 
The solid black lines are the upper and lower bounds ob- 
tained using the central values of the inputs, while the 
blue dashed and the green dashed-dot lines denote the 
weaker bounds obtained by including the errors using the 
first and the second method, respectively. 

In the rest of the analysis, the enlarged allowed inter- 
vals for the modulus due to the errors were obtained with 
the second method described above, which is actually the 
only one that allows for a proper inclusion of the correla- 
tions between the input values. 

We have devoted special attention to the input phase 
6(t) below ti n , as it plays an important role in the for- 
malism. We have used the phase-shift 5}(t) and the un- 
certainty from the CDF solution quoted in , presented 
in Eq. (jTSJ) and below it. The uncertainty of this phase is 
very small, even tiny at small energies, approaching half 
of a degree only at the upper end of the energy interval of 
interest. As a consequence, the effect of the phase uncer- 
tainty on the bounds is quite small. Alternatively, we have 
estimated the error of the phase as the difference between 
the parametrization (fT!?]) and the Roy solution recently ob- 
tained in [T]| . The results of this analysis are presented in 
Fig. [31 which shows the bounds obtained with the central 
input and the weaker bounds obtained by varying all the 
input quantities, using for the phase one of the two meth- 
ods of estimating the uncertainty discussed above. As can 
be seen, the two methods give practically the same results. 

We present in what follows a comparison of the bounds 
with the most recent experimental data from e + e~ anni- 
hilation [l|, 0, i, i, and r decays @. For the compari- 
son with e + e~ data we included in the calculation of the 
bounds the isospin breaking via the correction (|2T))) . 

The results shown in Figs. 0]- [HI are obtained with the 
weight p(t) = 1/t. In all of these, the central grey bands 
are defined by the upper and the lower bounds obtained 
with the central values of input quantities, while the ex- 
tended cyan bands are obtained by varying all the inputs 
within the errors as explained above. For clarity, we have 
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shown the bounds over the entire region of interest, and 
also separately in the low-energy and the resonance re- 
gions. 

As shown in Figs. [4] and H2 above the p peak the Belle 
data |(| lie within the grey band, therefore they are consis- 
tent with the bounds obtained with the central input. On 
the other hand, below the p peak the data are at the upper 
edge of the central band. The same features are observed 
in Figs. [7] and [SJ where we compare the bounds obtained 
with isospin breaking with the BaBaR [J, KLOE @, 3 
and CMD-2 @, || data. 

The behaviour below the p peak deserves attention, be- 
cause it could either indicate a problem in the data on the 
modulus, or suggest that our central input is not optimal. 
We first note that the bounds are quite sensitive to the 
value of (r 2 ) used as input. For illustration, in Figs. [TU1 and 
[IT] we show the bounds calculated with the central input 
used above, except for the choice (r 2 ) = 0.435 fm 2 instead 
of the central value quoted in (f^Tj) . This slightly higher 
value of the charge radius shifted upwards the bounds by a 
sizable amount in the resonance region. As a consequence, 
the data below the p peak are inside the allowed central 
grey band, but the data on the right side of the reso- 
nance fall now near the lower edge of the allowed region. 
Therefore, a very good consistency along the whole energy 
region is not achieved. With an even higher central charge 
radius, like (r 2 ) = 0.44 fm 2 , the lower bound at the left 
side of the resonance is situated above the data, so larger 
values of the charge radius lead to an even larger incon- 
sistency. Thus, the variation allowed for the charge radius 
is quite limited and cannot fully resolve the slight tension 
between the recent data on the modulus and the allowed 
inner band below the p peak. We must, however, empha- 
size that the data are fully inside the enlarged bands, and 
as a result there are no inconsistencies if we take into ac- 
count the errors on the input. 

One may also question the phase used as input. Indeed, 
while the F-wavc phase shift is obtained in @,H, H, EU EH 
by solving Roy equations for the nir amplitude, some input 
to these equations is borrowed in fact from the knowledge 
of the pion form factor itself. For instance, the value of the 
P-wave phase shift at the matching point of 0.8 GeV was 
taken in 0, l, E2 from a Gounaris-Sakurai parametriza- 
tion of the CLEO data on the form factor [28 1. A compar- 
ison of this paramctrization with the recent data shows 
that the latter are systematically larger than the older 
CLEO data (this is true for both the Belle data from r de- 
cay [|, and for the data [H, H 1, H from e+e" tt+tt", 
after applying the isospin correction). This implies that 
the input at the matching point 0.8 GeV in the Roy equa- 
tions might be improved with the new data on the form 
factor. 

In order to test the sensitivity of the bounds to the 
variation of the phase, we have repeated the calculation 
with the central phase from 0,113 increased (or decreased) 
by its quoted error (which is, in particular, of 2 deg at 
0.8 GeV). The results show that a higher phase leads to 
bounds shifted upwards, but it turns out that the change 
in the calculated bounds is not sufficient to improve the 



consistency, within the inner band, on the left side of the 
p peak. However, it appears that the bounds are more 
sensitive to the detailed shape of the input phase than to 
its overall magnitude. This aspect deserves further study. 

As we remarked above, for cross-checking the consis- 
tency of the data we must take into account the enlarged 
band obtained by allowing the inputs to vary within their 
errors. We have displayed the two kinds of bands since 
they impressively demonstrate how an improved input 
could lead to improved bounds. Furthermore, the inner 
band, obtained with the central values of the input, al- 
lows one to check whether a very good consistency holds 
(close central values are considered in general to indicate 
a better agreement of two determinations than a mere 
nonzero overlap of the intervals allowed by their error 
bars). Strictly speaking therefore, data lying at the edge 
of the inner band do not point at any tension between 
the data sets, and our results show in fact full consistency 
between the bounds and the new precise data on the mod- 
ulus. 

Turning to the low energy region, where the experi- 
mental errors are large, Figs.[5]and[8]show that the bounds 
are very stringent, both in the isospin limit and with 
isospin breaking, offering a more precise description of the 
modulus than experiment. This is perhaps the most inter- 
esting result of our analysis. Also, some tension between 
several points, especially from BaBar, and the bounds may 
be noted at low energies, despite the rather large experi- 
mental errors. 




Fig. 6. Bounds around the resonance region obtained with the 
weight p(t) — 1/t in the isospin limit. The bands are as in Fig. 
□ 



We have tested also other weights p(t) of the type 
given in ([3]). Some choices slightly improve the bounds in 
some energy regions, leading however to weaker bounds in 
other regions. Also, while one weight improves the upper 
bounds, other weights may lead to better lower bounds. 
In general, if we limit ourselves to bounds with a suitable 
decrease at large t such as to ensure a reliable calculation 
of the integral the differences in the resulting bounds 
are not very large. For illustration, we present below the 
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Fig. 7. Bounds obtained with the weight p{t) = 1/t with Fig- g Bounds around the resonance region obtained with the 



isospin breaking correction. The bands are as in Fig. [4] Also 
shown are the BaBar Q, KLOE 0,0] and the CMD-2 0, Q 
data. 



weight p(t) — 1/t with isospin breaking correction. The bands 
are as in Fig. [J] 



results for the weight p(t) = t x ' 2 /(t + 3), which leads to 
the value of the integral © given in Table [T] 
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Fig. 8. Low-energy profile of the bounds obtained with the 
weight p(t) = 1/t with isospin breaking correction. The bands 
are as in Fig. [4] 

In Figs. [T2] - HH we show the bounds obtained in the 
isospin limit, while Figs. [15]- [17] are obtained by including 
the isospin correction (|20[) . The central yellow bands are 
the allowed ranges for the modulus squared obtained with 
the central values of the input quantities, while the orange 
bands are obtained by varying all the inputs within the 
errors. For central values of the input this weight leads 
to upper bounds slightly lower than those obtained with 
p(t) = 1/t. As a consequence, the experimental data on 
the modulus lie now slightly above the central yellow band, 
both in the isospin limit and in the symmetry breaking 
case. Thus, the slight inconsistency remarked above for 
the bounds obtained with the weight pit) = 1/t is even 
more pronounced when the bounds are calculated with 
the weight p(t) = t 1 / 2 /{t + 3). However, the upper orange 
band is now wider and extends to higher values, so this 
weight finally gives weaker upper bounds than the prior 
one if the uncertainties are taken into account. On the 
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Fie. 10. Bounds on \F(t)\ 2 in the isospin limit, obtained with 
the central values on the input described in the text, except 
for (r 2 -) = 0.435 fm 2 instead of the central value in (|21|l . Also 
shown are the Belle data [jj. 



other hand, the lower bounds are slightly better in the 
resonance region and at higher energies than those given 
by the weight p(t) = 1/t, even when the uncertainties are 
taken into account. In the low-energy region, as for the 
previous weight, the bounds are very tight. 

Based on the above discussion, one can define "opti- 
mal" bounds by combining the results provided by differ- 
ent weights. For instance, we can take the upper bound 
given by the weight p(t) = 1/t, and the lower bound given 
by p(t) = t 1 ! 2 /{t + 3). These slightly improved bounds 
are presented in Figs. [18] - [23] both in the isospin limit 
and with isospin breaking corrections. We have performed 
the exercise mainly for purposes of illustration, as the im- 
provement is actually limited. Therefore, the comparison 
with the experimental data [II, LI, El III IE El leads to conclu- 
sions similar to those already formulated for the specific 
weights: above the p peak the data lie comfortably within 
the central band, below the p peak they are near the upper 
edge of the central band, while in the low-energy region 
the bounds are very strong, superseding the experimental 
data. 
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Fig. 12. Bounds on \F(t)\ 2 below t ln obtained with the weight 
p{t) = t 1/2 /(t + 3) in the isospin limit. The central yellow band 
is defined by the upper and the lower bounds obtained with 
the central values of the input quantities, the extended orange 
bands are obtained by varying all the inputs within the errors. 
Also shown are the Belle data [(J. 

5 Discussion and Conclusions 

Analyticity and unitarity have been used in a large num- 
ber of analyses of the pion electromagnetic form factor 
(for a list of references see t 2Ct] ) . A new formalism based 
on analyticity and unitarity, which exploits in an optimal 
way the precise information available on both the space- 
like and the timelike axes, was applied in fl9j and [20 | for 
constraining the shape of the form factor at t = and 
on the spacelike axis, respectively. In the present work we 
applied the same formalism for deriving upper and lower 
bounds on the modulus of the pion form factor on the 
unitarity cut below the lott threshold. 

The phase and the modulus are related by analytic- 
ity, but a complete reconstruction of the modulus requires 
the exact knowledge of the phase along the entire unitarity 
cut and the position of the possible zeros in the complex 
plane. As this information is not fully available, the cor- 
relation between the phase and the modulus is commonly 
investigated by means of specific parametrizations, or by 
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Fig. 14. Bounds around the resonance region obtained with 
the weight p(t) = t 1 ^ 2 /(t + 3) in the isospin limit. The bands 
are as in Fig. 1121 



making ad-hoc assumptions in dispersive representations 
of the Omnes type. 

Our aim was to check the consistency between the 
phase and the modulus of the pion form factor in a model 
independent framework. We have used as input the phase 
in the elastic region, known via the Fermi- Watson theo- 
rem from the phase shift of the P-wave of 7T7r scattering. 
The lack of information on the phase above the inelastic 
threshold and the zeros in the complex plane has been 
compensated by some experimental information and con- 
servative assumptions on the modulus above this thresh- 
old. We have used the BaBar data on the modulus [l[ in 
the range between the lot: threshold and 3 GeV, and very 
conservative assumptions above 3 GeV. Using modulus in- 
formation above ti n and the phase below i; n , it is possible 
to calculate upper and lower bounds on the modulus in 
the elastic region. 

In order to reduce the influence of the intermediate en- 
ergies, where the modulus is not known, we implemented 
the information on \F(t) \ above t- m through the L 2 integral 
condition $2$, instead of imposing this knowledge point- 
wise, at each t. Then the bounds are inevitably weaker, 
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but for suitable choices of the weight pit) they are not sen- 
sitive to the assumptions on the modulus at intermediate 
and high energies. We also included information on the 
normalization and charge radius at t = and the values 
at some spacelike momenta, available from recent experi- 
ments [ntm. 

Wc mention that, for a given input, the bounds are op- 
timal and imply no loss of information as they are obtained 
by solving exactly a functional interpolation problem. As 
seen from the formalism presented in Sec. [2j the results 
do not rely on specific paramctrizations. Furthermore, the 
knowledge of the phase of the form factor in the inelastic 
region is not required. Indeed, although some intermedi- 
ate quantities entering the solution of the problem involve 
the phase chosen as an arbitrary (smooth) function above 
ti n , it was proven and checked numerically [l7| that the 
dependence is exactly compensated among the relevant 
factors. 

Wc have compared the bounds obtained from our anal- 
ysis with the precise measurements of the modulus in the 
elastic region from recent experiments on e + e~ annihila- 
tion 

0, H,ll 0, II and hadronic r decays @. Of course, we 
could use as input also one of the measured values on the 
modulus below the inelastic threshold. This would lead to 
a considerable strengthening of the bounds. However, as 
the four sets of data are not perfectly consistent among 
them, the choice of the input from one set would intro- 
duce a certain bias, which we prefer to avoid in this first 
analysis. Since the problem is of great interest, it will be 
studied in the near future. 

The bounds obtained from our analysis in the isospin 
symmetry limit can be compared directly with the recent 
experimental data on the modulus available from hadronic 
r decays Q. For the comparison with the data obtained 
from e+e" annihilation 0,H i, 0, Q , isospin breaking was 
included in the formalism through the phase modification 
(j2TJ)) . In both cases, we have obtained a very good consis- 
tency above the p peak, where the data lie well within the 
inner bands derived with the central values of our input. 

On the other hand, below the p peak the data are 
situated at the upper edge of the central region. As this 
may indicate that our input is not optimal, we briefly an- 
alyzed in Sec. [4] the possible changes in the input that can 
improve the consistency. The results show that it is not 
possible to achieve a very good consistency at all energies 
by tuning the value of the charge radius: with a slightly 
larger input value, (r 2 ) = 0.435 fm 2 instead of the central 
value quoted in (|2ip . wc obtain good consistency below 
and around the p peak, but the agreement is worsened at 
higher energies. Furthermore, higher values of the charge 
radius are even less acceptable, as they lead to central 
bands in disagreement with the experimental data. 




1.3 0.4 0.5 0.6 0.7 0.8 0.9 



t" 2 [GeV] 

Fig. 15. Bounds obtained with the weight p(t) = t 1/2 /(t + 3) 
with isospin breaking correction. The bands are as in Fig. Q2] 
Also shown are the BaBar [J], KLOE 0, 0| and the CMD-2 
H H data. 
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Fig. 16. Low-energy profile of the bounds obtained with the 
weight p{t) — t 1//2 /(i + 3) with isospin breaking correction. The 
bands are as in Fig. 1121 



Concerning the input phase, preliminary studies show 
that the bounds are quite sensitive to the detailed shape 
of the phase in the elastic region. We recall that the in- 
put phase was obtained via the Fermi- Watson theorem 
from the P-wavc phase shift of tttt scattering. This phase 
shift was determined by solving Roy equations for 7T7T am- 
plitude, using however as input the phase at a matching 
point from phcnomcnological representations of the form 
factor. Thephase used in the present analysis is taken 
from 0, 0, ©, 03 (HI, which used previous measurements 
of the modulus (in particular, Refs. @, H, LUt took the in- 
put value of the P-wave phase shift at the matching point 
0.8 GeV from a parametrization of CLEO data [28]). Of 
course, one can adjust the input value until the phase ob- 
tained by solving Roy equations leads to upper and lower 
bounds consistent with the more precise recent data on 
the modulus. This more elaborate analysis is also a project 
for future work. Taken together with other recent investi- 
gations [23|, [24| , such a study will contribute to a better 
understanding of the form factor in the clastic region, cru- 
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Fig. 17. Bounds around the resonance region obtained with 
the weight pit) = t 1 ^ 2 /(t + 3) with isospin breaking correction. 
The bands are as in Fig. 1121 
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Fig. 18. Optimal bounds on |F(t)| 2 below t- m in the isospin 
limit, obtained with the weight p(t) = 1/t for the upper bounds 
and pit) = t 1 ^ 2 1 it + 3) for the lower bounds. The central grey 
band is defined by the upper and the lower bounds obtained 
with the central values of the input quantities, the extended 
violet bands are obtained by varying all the inputs within the 
errors. Also shown are the Belle data 0]. 

cial for the precise evaluation of the hadronic contribution 
to the anomalous magnetic moment of the muon [2^, [26| . 

The considerations above, based only on the inner band, 
illustrate mainly the sensitivity of the results to the vari- 
ation of the input. For checking the consistency we must 
of course take into account the enlarged bands, obtained 
by allowing the input to vary within errors. From the re- 
sults presented in Section [4] we conclude that there are no 
inconsistencies between the bounds and the recent exper- 
imental data on the modulus above 0.5 GeV. 

At low energies the bounds derived in the present work 
are very stringent, leading to allowed ranges for the mod- 
ulus smaller than the quoted errors of the experimental 
measurements. This is one of the most important result of 
our work, showing that the knowledge of the form factor at 
low energies, specifically below 0.5 GeV, can be improved 
by exploiting in an optimal way the phase in the elastic re- 
gion and conservative information on the modulus above 
the inelastic threshold. Also, below 0.4 GeV the results 
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Fig. 19. Low-energy profile of the optimal bounds in the 
isospin limit. For details see caption to Fig. 1181 
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Fig. 20. Optimal bounds around the resonance region in the 
isospin limit. The bands are as in Fig. [18] 
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Fig. 21. Optimal bounds with isospin breaking correction. The 
bands are as in caption to Fig. 1181 Also shown are the BaBar 
Q, KLOE and the CMD-20,01 data. 

reveal some inconsistencies between the bounds and the 
experimental data on the modulus, especially from BaBar, 
despite the large errors on these data. 

We recall that the bounds are obtained with a defi- 
nite range for the charge radius, given in (|2ip . and are 
sensitive to this input. As the adopted range (|2"Tj) leads to 
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Fig. 22. Low-energy profile of the optimal bounds with isospin 
breaking correction. The bands are as in Fig. 1181 
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Fig. 23. Optimal bounds with isospin breaking correction 
around the resonance region. The bands are as in Fig. 1181 

consistency for data on modulus at higher energies, but to 
some inconsistencies at low energies, the results indirectly 
signal some inconsistencies between data themselves. As 
mentioned in the Introduction, the formalism presented in 
this paper can be exploited also the other way round, for 
finding bounds on the pion charge radius from the precise 
measurements on the modulus below the uitt threshold. 
This analysis will be performed in a future work. 

We finally mention that further improvements are ex- 
pected if very precise determinations of the form factor 
from ChPT, lattice simulations or experiment are made 
available. Indeed, as we have shown, precisely known in- 
put values of the form factor inside the analyticity domain 
lead to a considerable strengthening of the bounds. There- 
fore, the present formalism is a useful framework, open to 
further developments, for analyticity tests and precise pre- 
dictions on the pion form factor at low and intermediate 
energies. 
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